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We present results of the detailed Monte Carlo calculation of the rates of double-bang events 
in 1 km'^ underwater neutrino telescope with taking into account the effects of r-neutrino 
propagation through the Earth. As an input, the moderately optimistic theoretical predictions 
for diffuse neutrino spectra of AGN jets are used. 
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I. INTRODUCTION 

During these years there have been many discussions 
[1,2] on the possibihties of very high energy r-neutrino 
detection by large neutrino telescopes, which would prove 
oscillations of muon neutrinos of astrophysical origin 
("astrophysical long basehne experiment" [3]). 

It is known now that neutrino oscillations definitely 
exist. Nevertheless, the detection of astrophysical r- 
neutrinos is still very interesting and useful, e.g., for a 
better understanding of neutrino properties as well as 
properties of astrophysical sources [4] . 

Intrinsic flux of r-neutrinos from a typical source in 
which there is acceleration of protons and production of 
neutrinos through pp- and p7-reactions is very small. The 
dominant channel of r-neutrino birth is through the in- 
clusive production of charmed mesons [5] , 
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and the smallness of Vr-Q-\i^ follows from the relations 
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Therefore, the production of Vt at source is negligible, 
and if all muons can decay, it can be expected for the 3 
flavors the following proportion: 
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Assuming maximal atmospheric mixing and Uez = 0, 
that is: 



sin ^1 3 = , sin 26*; 
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one obtains, after averaging over many oscillations, the 
ideal equipartition between neutrino flavors in the astro- 
physical neutrino flux at Earth, 
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independently on the value of the solar mixing angle 9i 2 . 
Moreover, in the case when there is maximal atmospheric 
mixing and [/e 3 = 0, for any proportions between flavors 



at the source one has, after propagation, the equality 

•/-M = 't'r [4]. 

The experimental verification of the relation (f)^ — 
by detection of astrophysical diffuse fluxes of Vr and z/^ 
with flavor identification capabilities is very important 
because the inequality 0^ ^ (pr is predicted in some ex- 
otic scenarios, e.g. [4]: i) if CPT invariance is violated 
in the neutrino sector and the neutrino and antineutrino 
mixing matrices do not coincide, i/)^, in general, is not 
equal to ii) if neutrinos can decay [6] en route from 
the sources to the Earth, the neutrino fiuxes and flavor 
ratios are very sensitive to uncertanties in the neutrino 
mixing matrix and strongly depend on the hierarchy of 
neutrino masses. So, the inequality 0^ ^ (pr would reveal 
unconventional neutrino physics. 

Inside of some astrophysical and cosmological objects 
there are processes in which very energetic quarks and 
gluons fragmenting into jets of hadrons are expected to 
be produced. Even the usual inclusive spectrum of high 
energy pp- or p7-collisions always contains the jet com- 
ponent (described well by perturbative QCD). Jets can 
appear in decays of (hypothetical) supermassive parti- 
cles which, in turn, are produced in decays of topolog- 
ical defects [7], in processes of superheavy dark matter 
annihilation inside of stars [8], in a process of Hawking 
evaporation of primordial black holes [9] . It was shown in 
[10] that decays of top quarks contained in the jets lead 
to production of r-neutrino flux from the jets, which is of 
the same order as Vp,, Ve fluxes (at large neutrino energies, 
close to the value of the jet mass). This means that the 
intrinsic r-neutrino flux from the sources in which the jet 
phenomena are essential is, in general, not suppressed (in 
comparison with intrinsic v^, v^, fluxes). 

We studied in this paper the possibility of detection of 
extragalactic diffuse r-neutrino fluxes by large (~ Ikm'^) 
underwater neutrino telescopes. In Sec. II we discuss the 
general problems connected with the detection of very 
high energy r-neutrinos {E^ > lO^GeV). In Sec. Ill we 
present the diffuse neutrino spectra for the AGN jets cho- 
sen for the numerical calculations. Details of the calcu- 
lations, results and conclusions are given in Sec. IV. 
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II. DETECTION OF r-NEUTRINOS 

The main feature of charge current interactions of r- 
neutrinos is the fact that r-leptons decay long before they 
lose a large fraction of their energy. This leads to the 
absence of absorption of r-neutrinos during their propa- 
gation through the Earth and , as a consequence, to the 
characteristic modification of the neutrino spectrum after 
crossing the Earth (there is a pile-up of events with ener- 
gies around ~ lOOTeV, if the incident neutrino spectrum 
is not too steep). At neutrino energies ^ 10^ — lO^GeV 
T-lepton track-length becomes larger than ~ 100m (in 
water), it becomes possible to identify r-neutrino in- 
teractions in large neutrino telescopes by selection of 
"double-bang" or "lollipop" events [1]. A double-bang 
event consists of two showers and a r-lepton track con- 
necting them. The first (hadronic) shower is initiated by 
the charge current interaction of r-neutrino and the sec- 
ond one (hadronic or electromagnetic) is produced by the 
T-lepton decay. In lollipop events one of these showers 
is not detected (e.g., if the r-neutrino interaction takes 
place outside of the detector) . 

The detection of events with double-bang structure 
would reliably indicate the appearance of r-neutrinos in 
the detector because this signature is unique. However, 
the probability of observing a double-bang event Pdb 
decreases with neutrino energy for > lO^GeV (due 
to r-lepton range increase). The lollipop structure is not 
so spectacular although the observation probability Pl 
of lollipop events is at E^, ^ lO^GeV considerably (by a 
factor of 10 [11]) larger than Pdb- 

Diffuse astrophysical and cosmological neutrino fluxes 
in the ultra-high energy region, E^, > lO^GeV, (GZK 
neutrinos (see, e.g., [12]), hypothetical neutrino fluxes 
from topological defects [13], from Z-bursts [14] etc.) are 
too small and inaccessible for a study by underwater neu- 
trino telescopes. Registration of r-neutrinos of such en- 
ergies can be accomplished by detecting the air show- 
ers caused by r-leptons produced by neutrino-nucleon in- 
teractions far away from the detector ("Earth-skimming 
idea" [15]). The huge modern air shower detectors (e.g., 
the HiRes [16] or the Pierre Auger detector [17]) could be 
suitable for such experiments, but expected event rates 
are low. 



III. DIFFUSE NEUTRINO SPECTRA FROM 
AGN JETS 

The concrete calculations were carried out with two 
theoretical diffuse neutrino spectra from AGNs. 

1. The curve denoted by "MPR(limit)" in Fig.l is the 
generic upper bound on the diffuse neutrino spectrum 
from AGN jets obtained by Mannheim, Protheroe, and 
Rachcn [18]. This limit was found using the assumption 
that the AGN source is not completely optically thin 
for a cosmic ray flux: due to opacity effects there is a 
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FIG. 1. Different theoretical predictions for diffuse neu- 
trino (j/p+i)^) spectra from AGNs. MPR(limit): upper bound 
for flux from AGN jets [18], P: flux for proton blazar 
model of [21], SSQC: prediction of flux from radio-quiet 
AGNs [22]. 

spectral break between 10^ GeV and 10^^ GeV in the 
escaping cosmic ray spectrum (the resulting cosmic ray 
spectrum is a combination of several model spectra with 
different values of the spectral break and with the fixed 
Emax {Emax — 10^^ GcV), and the normalization of 
the resulting spectrum is such that the cosmic ray inten- 
sity does not exceed the proton spectrum estimated from 
observations). Neutrinos are mostly produced in decays 
of mesons from accelerated proton interactions with syn- 
chrotron photons emitted by high energy electrons in the 
jets. It is noted in [18] that various models of AGN jets 
(see, e.g., [19-21]) predict diffuse neutrino spectra which 
are compatible with the MPR limit (at the important 
neutrino energy range below 10'' GeV). 

2. The model by Protheroe [21] (the curve denoted by 
"P" in Fig.l) is the optically thick proton blazar model, 
in which it is assumed that the target for the pion pro- 
duction is provided by the thermal UV photons emitted 
from the accretion disc rather than by the synchrotron 
photons produced in the jets. The diffuse neutrino spec- 
trum is normalized to the experimental data on 7-ray 
background. It is seen from the figure that the predic- 
tion of this model is, indeed, not in the contradiction with 
the MPR bound, if E^ < lO'^ GeV. The proton blazar 
model, developed by Mannheim [19] (in which acceler- 
ated protons interact with the synchrotron radiation of 
jet electrons) predicts very similar diffuse neutrino spec- 
trum. 

One should note that the prediction given by Steckcr 
and Salamon [22] for diffuse neutrino spectrum from 
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radio-quiet AGNs is much more optimistic (SSQC curve 
in Fig.l) at neutrino energies E^, ^ 10^ — 10^ GeV. Their 
model, however, predicts a non-thermal spectrum of X- 
rays at energies below ~ 500 KeV{these X-rays result 
from reprocessing of gamma-rays from pion decays by 
electromagnetic cascades due to the high photon density 
near the black hole) while the observed X-ray spectra of 
radio-quiet AGNs (as well as the spectrum of the diffuse 
X-ray background) turn over steeply below ^ 100 KeV 
with no sign for a nonthermal component [23]. Therefore, 
the existing X-ray data cannot be used for a normaliza- 
tion of the radio-quiet AGN neutrino spectrum (authors 
of [22] normalized their spectrum, rather arbitrarily, on 
30% of the diffuse X-ray background). This is the rea- 
son why we did not use SSQC flux in our calculations. 
Besides, SSQC neutrino flux is now almost excluded by 
measurements on AMANDA [24] and Baikal [25] neutrino 
telescopes. 

IV. RESULTS AND DISCUSSIONS 

First estimates of double-bang event rates in km'^ neu- 
trino telescope were done in [26], for down-going r- 
neutrinos only. Detailed calculations of r-neutrino prop- 
agation through the Earth carried out by Dutta ct al. [27] 
(using power law incident neutrino spectra) showed that, 
in spite of the absence of neutrino absorption, the contri- 
bution of lower hemisphere events to the total event rate 
in a telescope cannot be too large. This is due to a strong 
decrease of neutrino energies during propagation. It was 
concluded in [27] that the most promising way to detect 
r-neutrino astrophysical fluxes is an analysis of shower 
events in the telescope. 

In this work we present results on Monte Carlo calcula- 
tions of rates of double-bang events in a 1 km^ telescope 
for both hemispheres using the diffuse neutrino spectra 
shown on Fig.l. We simulate both DIS CC and NC 
neutrino interactions for all flavors in the energy range 
lO^GeV < < lO'^GcV using the CTEQ3_ DIS struc- 
ture functions. For r-lcpton propagation we used the 
MUM code [28] that was updated accounting for the cor- 
rections for photonuclear interactions of r-leptons [29]. 
Some details of the calculations and the resulting neu- 
trino fluxes of all flavors after propagation through the 
Earth were presented in our previous work [30] . 

The rate of totally contained double-bang events is 
given by the formula 

/I />oo 
/ V,ffiE,^,0)I{E,^,0)x 
-1 JEmin 

a^^{E^^)dE^^d{cose). 

Here, Na is the Avogadro number, p is the density of 
the medium (we use p = Ig cm~^ which is close to sea 
water /ice density), I{Ev^,6) is the diflFerential Ur flux 



entering the detector at zenith angle 9, is the to- 

tal deep inelastic neutrino cross section (charge current 
part). K:// is the eff'ective volume given by the relation: 

K// (£;.., e) - S^{Q){L - Rr{E,^)), 

where Sp is the projected area for tracks generated 
isotropically in azimuth at the flxed 9 directions on a 
detector of parallelepiped form, L is the geometrical dis- 
tance between entry and exit point, Rr is the r lep- 
ton range. The minimum value of the neutrino energy, 
Emin = 2 • 10^GeV,corresponds to a r-lepton range of 

The double-bang topology can be considered as a back- 
ground free signature of r-neutrino appearance [1] since 
there should be no atmospheric muons which could pro- 
duce in the detector 2 showers with comparable amounts 
of photons. Nevertheless, this reasonable assumption 
needs to be verified through a full simulation. 

Concrete calculations were done for two geometries of 
the neutrino detector: 

1x1x1 km^ (IceCube-like [31]), 
1.4 X 1.4 X 0.6 km^ (NEMO-like [32]). 

In Table 1 the corresponding double-bang rates 
for lower, upper and both hemispheres are shown. 



TABLE 1. Number of totally contained double- 

bang events in km^ detector per year. 



Spectrum 


Rate {N-2n/N2n/Ni^) 




NEMO-like 


[21] 


1.0 2.1 3.1 


[18] 


1.4 3.1 4.5 




IceCube-like 


[21] 


0.7 1.6 2.3 


[18] 


1.0 2.3 3.3 



Calculated values for upper hemisphere are 3 — 6 times 
lower compared to [26] since more optimistic predictions 
for diffuse neutrino fluxes were used there. 

Finally, we found using the theoretical predictions in 
[18,21] for AGN diffuse neutrino fluxes that in a 1 km^ 
neutrino telescope one can expect a marginally observ- 
able rate of double-bang events from r-neutrinos. 
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